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I.  INTRODUCTION 


A  substantial  amount  of  research  has  been  conducted  on  the  stability  of 
the  motion  of  spinning,  liquid-filled  containers.  Areas  of  interest  can  be 
quite  diverse.  Two  such  examples  are:  the  geophysical  problem  of  the  motion 
of  the  fluid  core  of  the  earth  or  the  flight  stability  of  a  spin-stabilized 
liquid  payload  artillery  projectile.  Most  experience  in  rotating  liquid  flows 
is  for  cases  where  the  rotational  force  is  much  larger  than  the  viscous  force, 
i.e.,  high  Reynolds  number  (Re)  flows.*  This  paper  deals  with  a  series  of 
experiments  that  were  conducted  for  5  <  Re  <  12,000.  This  range  of  Re  has 
previously  not  been  examined. 

Stewartscn  treated  the  stability  of  the  motion  of  a  liquid-filled  top.1 
In  this  analysis,  the  inviscid  liquid  was  assumed  to  be  spinning  as  a  quasi- 
rlgld  body  with  the  angular  velocity  of  the  top.  It  was  determined  that  the 
yawing  motion  of  the  top  could  excite  the  natural  frequencies  of  oscillation 
of  the  liquid  (the  eigenf requencies)  and  a  moment  due  to  the  pressure  distur¬ 
bances  within  the  liquid  could  destabilize  the  motion  of  the  top.  Wedemeyer2 
incorporated  viscosity  into  the  Stewartson  theory  through  the  use  of  boundary 
l^yer  type  corrections,  and  the  Stewartson-Wedemeyer  theory  has  compared 
favorably  with  gyroscope  experiments  conducted  by  Karpov.3  Recent  analyses  by 
Murphy4  and  by  Gerber,  et  al,  5  represent  more  complete  and  rigorous  models 
than  the  Stewartson-Wedemeyer  theory,  but  all  of  these  treatments  are  still 
limited  to  "high  Reynolds  numbers"  due  to  viscous  corrections  at  the  walls  of 
the  container. 


*  Reynolds  number  is  defined  as  afyv  (whore  a  is  the  radius  of  a  oylindrioal 

♦ 

aavity,  $  is  t ho  spin,  and  v  is  the  kinematic  viscosity  of  the  liquid )  and 
is  physically  the  mtio  of  the  rotational  force  to  the  viscous  force .  Thie 
definition  is  simply  the  inverse  of  the  Kkmn  number  which  is  typically 
used  in  analyses  of  rotating  liquid  flows. 

1.  K,  Stewartson,  “On  the  Stability  of  a  Spinning  Top  Containing  Liquid,  " 
Journal  of  Fluid  Mechanics,  Vol .  S,  1959,  pp.  577-593. 

2 .  R.  H.  Wedemeyer,  “Viscous  Correctiom  to  Stewartson' a  Stability  Criteri¬ 
on,  "  8RL  Report  No.  1325,  June  196 6  (AD  489687), 

3.  8.  G.  Karpov,  “Liquid- Filled  Gyroscope:  The  Effect  of  Reynolds  Number  on 
Resonance,  “  BRL  Report  No.  1302,  October  1965  (AD  479430). 

4.  Charles  H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  With  a  Viacoue 
Liquid  Payload,  “  BRL  Memorandum  Report  No.  ARBRL-MR-03194,  August  1982 

(AD  A1286W. 

5.  Nathan  Gerber  and  Raymond  Sedney,  “Moment  on  a  Liquid-Filled  Spinning  and 
Nutating  Projectile:  Solid  Body  Rotation,  “  BRL  Technical  Report  No. 
ARBPL-TR-024 70  (AD  A125332).  ' 
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for  low  Reynolds  number  rotating  flows,  D'Amico  and  Miller  observed  large 
yaw  and  despin  moments.6  Miller  conducted  laboratory  experiments  using  a 
fixed  precession  angle  spin  fixture  and  determined  despin  moments  over  a  wide 
range  of  Reynolds  number  (Figure  l).7  The  maximum  despin  moment  is  at  Re  = 

10.  D'Amico  fired  yawsonde-instrumented  projectiles  which  produced  rapid  yaw 
growth  and  violent  despin  histories  for  approximately  that  same  Reynolds 
number.8*9  This  report  describes  a  series  of  laboratory  experiments  where  the 
1 i quid -induced  yaw  moment  was  measured  as  a  function  of  the  yawing  frequency 
of  a  gyroscope.  The  spin  of  the  gyroscope  was  held  constant  (no  despin  mo¬ 
ments  will  be  reported),  and  the  yawing  motion  of  the  gyroscope  was  less  than 
5  degrees.  The  Reynolds  number  range  for  these  data  are  5  <  Re  <  12,000. 
Comparisons  of  liquid-induced  yaw  moment  coefficients  for  flight  and  gyroscope 
tests  are  made.  Correlations  for  the  liquid-induced  yaw  moment  with  perti¬ 
nent  dimensionless  groups  are  made.  Also,  a  lower  bounds  in  Reynolds  number 
was  tentatively  established  for  the  application  of  viscous  corrected  theories 
such  as  Murphy4  and  Gerber,  et  al.5 


II.  LABORATORY  GYROSCOPE 

A.  Simulation  of  a  Spin-Stabilized  Projectile  With  a  Gyroscope 

A  freely  giinballed  gyroscope  similar  to  that  used  by  Karpov3  was  used  to 
measured  the  yaw  moments  produced  by  a  highly  viscous,  rotating  liquid.  A 
short  explanation  of  the  dynamical  behavior  of  the  gyroscope  will  aid  in  the 
interpretation  of  the  test  data  and  results.  A  spin-stabilized  projectile 
experiences  a  complicated  angular  motion  during  free-flight.  Only  certain 
aspects  of  this  free-flight  motion  will  be  simulated  by  the  gyroscope.  Murphy 
has  analyzed  the  motion  of  a  spin-stabilized  projectile10  For  the  case  of  a 
projectile  which  does  not  contain  a  liquid  and  where  the  effects  of  drag  are 
neglected,  the  angular  motion  is  ascribed  by 


6.  V.  P.  D'Amico  and  M.  C.  Miller,  " Plight  Instability  Produced  by  a  Rapidly 
Spinning,  Highly  Viscous  Liquid,  11  Journal  of  Smoaora,ft  and  Rockets.  Vol . 
16,  No.  1,  January-February  1979,  pp.  62-64 . 

7,  Miles  C.  Miller,  " Flight  Instabilities)  of  Spinning  Projectiles  Raving 
Non-Rigid  Payloads,  "  Journal  of  Guidance  Control  and  Dynamics,  Vol .  5, 
Maroh-Apnil  1962,  pp.  151-157, 

J.  V.  P.  D'Amico,  V.  H.  Clay,  and  A.  Mark,  "Diagnostic  Teats  for  Vick-Type 
Payloads  and  High  Viscosity  liquids,  "  BRL  Memorandum  Report  No.  ARBRL-MR- 
02913,  April  1979  (AD  A072812). 

9.  V.  p,  D'Amico  and  V.  H.  Clay,  * High  Viscosity  Liquid  Payload  Xamonde 
Data  for  Small  Launch  Yaws,  "  SRL  Memorandum  Report  No.  ARBRL-MR-03929, 
June  1980  (AD  A08F411). 

10.  Carles  H.  Murphy,  "Free  Flight  Motion  of  Syrmstrio  Missiles,  "  BRL  Report 
No.  1216,  July  1963  (AD  442757). 
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(1) 


5  +  (H  -  i  a4>)  5  -  (M  +  io$T)  5=0 


where  5  is  the  complex  yaw  within  a  nonrolling  coordinate  system  and 


H  =  damping  moment  a  (pSi/2m) 


cNa  -  ky'2 


(S +  Ch;) 


M  =  static  moment  =  (pSz3/2Iy)  CM  (V/t): 


f  =  Magnus  moment  =  (pS£/2m) 


+  kx  2  C| 


V 


t*/t) 


$  =  spin  rate 


0  =  IX/Iy 


If  an  epicyclic  motion  is  assumed  for  the  etipty  projectile,  then 


i  i$2 

5  *  K^e  ♦  K2e 


(2) 


where  ordinarily 

G-i  tit 

Kj  «  Kj  e  J  J  ,  j  *  1,  2 

h  “  ho *  Tj¥ 

Tj  =  (a/2)  (1  t  y]\  -  I/s“) 


(3) 

(4) 

(5) 
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where  is  the  gyroscopic  stability  factor,  Ix  and  Iy  are  the  axial  and 
transverse  moments  of  inertia,  and  x^  and  x2  are  the  fast  and  slow  yaw  fre¬ 
quencies.  (A  spin  stabilized  projectile  has  two  modes  of  precession  or  yaw.) 
Explicit  forms  for  the  yaw  growth  rate  (e.-)  and  the  gyroscopic  stability 
factor  are  J 


ej  =  "  (Tj^  "  0  T)/[(2tj  -  o)  ( 4>) 3  (6) 

sg  =  o2  4>2/ 4M  (7) 

A 

A  gyroscope  has  a  gravity  moment  (G)  that  acts  in  a  similar  fashion  as 

the  static  moment  M.  Experience  with  the  gyroscope  has  shown  that  both  pro¬ 
cessional  modes  will  damp  when  the  center  of  mass  is  slightly  below  the  gimbal 

A 

pivots.  For  this  configuration  G  <  0.  The  direction  of  ^  will  be  the  same 
as  $,  while  the  direction  of  will  be  opposite  to  $.  The  gyroscope  sta¬ 
bility  factor  will  be  negative  since  G  <  0.  However,  | j  »  1  and  for  the 
empty  gyroscope. 


a  o 


and 


t2 


3  0 


(8) 


Typically,  the  gyroscope  is  operated  such  that  0.Q4  <  <  0.10  and  | 

<  0.005,  which  corresponds  to  the  conditions  set  out  in  Equation  (8). 

A  form  for  the  liquid  moment  has  also  been  suggested  by  Murphy.1* 


■if  +  1  a  \  a2$2 


T1  CLH1  K1 


Hz 

6  +  T2  CLM2  K2  e 


(9) 


where  the  liquid  moment  was  scaled  by  the  liquid  mass  (mj,  cylinder  radius 

(a),  and  the  spin  {$).  This  internal  liquid  moment  can  be  incorporated  into 
the  projectile  description  provided  by  Eq.  (1), 


t  >  (H  -  io*$)  c  .  (h  +  io^  f)  l 


t  i2  ("H.a2/ly)  [  >i  CLHi  K, 


1  $1 

e  +  T2  CLH2  k2  e2 


i  $2  J 


(10) 
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;> Tcf^.^^r^g-Kiry^t  !>^?''-;-'r' ■  wt » i-  '■-'*■  --  ■>  -.■ ■  -. 


Several  assumptions  valid  for  the  gyroscope  can  be  made  to  simplify  this 
equation.  First,  damping  rates  for  the  empty  gyroscope  are  very  small  (typi- 

cally  s j  <  -  5  x  10"'3).  Hence,  the  only  moment  acting  on  the  gyroscope  will 
be  that  produced  by  the  liquid.  Second,  the  yawing  motion  will  be  assumed  to 

~  i  4>i 

consist  of  only  the  fast  yaw  mode;  i.e, s  £  =  K^e  .  Third,  CLM  will  ulti¬ 
mately  depend  upon  many  variables  (K^,  tj,  e^.  Re,  c/a),  but  a  simple  form 
from  Reference  4  will  be  used  when  the  spin  is  positive. 

clMj  s  cLSM.j  +  1  cLIMj  (n) 

CLSM  is  a  side  moment  coefficient  that  controls  the  projectile  system  yaw 
damping,  while  is  an  in-plane  moment  coefficient  that  principally 
modifies  the  frequency  of  motion. 

General  expressions  for  yaw  frequencies  and  yaw  growth  (or  damping)  rates 
can  be  obtained  from  Equations  (10)  and  (11): 


tj  “  un)  cfj  -  (-uJ  ^ f]  -  (i/sg)  ] 


(12) 


for 


fj  «  1  +  (mL  a2/Ix)  CLIM^ 


and 

ej  -  (\a2/lx)  (2tj/o  -1)  1  CLSHj 


(13) 


Equations  (12)  and  (13)  can  be  inverted  to  obtain  expressions  for  CLIM 

and  Glsm  based  upon  experimentally  determined  values  of  ej  and  xy  During 

j 


'j 


the  course  of  the  experiments  it  was  observed  that  t  for  the  liquid-filled 
gyroscope  was  within  1*.  of  c  for  the  empty  gyroscope.  Inis  small  difference 
is  the  effect  of  C|_IM  on  t.  Attempts  were  made  to  utilize  this  small  differ¬ 


ence  between  the  liquid-filled  and  empty  coning  frequencies  to  establish  CLIM 
values.  It  was  found  that  the  frequency  measurements  were  not  sufficiently 
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accurate  to  predict  Cum*  Equation  (13)  will  be  used  to  determine  C^SM. 

The  stability  of  the  gyroscope  Is  controlled  by  the  liquid,  and  the  motion  is 
dominated  by  the  fast  yaw  mode.  The  yawing  motion  is  essentially  a  single 
'  mode  coning  motion.  As  such,  the  dimensionless  fast  yawing  frequency  will  be 
referred  to  as  the  coning  frequency,  i.e.,  t  *  t^,  and  the  yaw  growth  rate  of 

the  fast  mode  Is  redefined  as  =  e.  For  t  *  o  and  H  »  f  »  0,  then.  Equation 

(13)  becomes 


CL$M  3  c/ (mi_a2/Ix)  • 


(14) 


The  above  formulation  also  assumes  that  c  Is  relatively  small.  For  the  exper¬ 
imental  data  presented  here,  e  <  0.02,  and  Equation  (14)  can  be  used. 

B.  Operation  of  a  Liquid-Filled  Gyroscope 

Cylindrical  containers  are  located  within  the  rotor  of  a  freely  glmballed 
gyroscope  (Fig.  2).  The  Inner  glmbal  supports  the  rotor  and  the  test  cylin¬ 
ders.  Weights  can  be  located  on  the  rotor  to  control  Ix.  Weights  can  also  be 

attached  to  the  top  or  bottom  of  the  Inner  glmbal.  During  the  course  of  a 
test  run,  a  DC  motor  (located  beneath  the  Inner  glmbal)  drives  the  rotor  at  a 
constant  rate.  The  positions  of  the  weights  located  on  the  Inner  glmbal  can 
be  moved  to  vary  Iy,  and  these  weights  are  used  to  control  Ix/Iy,  which  from 

Eq.  (8)  provides  control  of  the  coning  frequency.  The  center  of  gravity  of 
the  rotor/inner  glmbal  Is  selected  to  be  slightly  below  the  support  pivots. 

The  pivots  consist  of  crossed  spring  leaves  (so  called  flexural  pivots)  in¬ 
strumented  with  strain  gages.  With  the  aid  of  usual  bridge  circuit  tech¬ 
niques,  strain  Is  calibrated  as  a  function  of  deflection.  The  response  of  the 
measurement  system  is  linear  for  deflections  up  to  ten  degrees. 

The  yaw  of  tha  gyroscope  Is  the  angle  between  the  vertical  and  the  spin 
axis  of  the  rotor.  The  liquid  Is  allowed  ample  time  to  achieve  a  state  of 
rigid  body  rotation  prior  to  the  rotor  being  released.  The  Inner  glmbal  Is 
released  from  the  vertical  position  without  aqy  initial  disturbance.  Unstable 
motion,  If  it  occurs.  Is  normally  self  starting  from  the  vertical  position. 

The  yawing  histories  for  each  set  of  experiments  are  digitized  and  passed  to  a 
computer.  The  computer  processes  the  yaw  record  to  determine  a  Kj  history. 

An  initial  amplitude  Is  determined,  K10,  and  Is  used  as  a  reference  by  which 

all  successive  Kj  values  are  scaled.  The  ratio  [In  (Kj/K10)]  Is  displayed 

versus  the  time.  A  sample  reduction  for  In  (Kj/Kjq)  versus  time  Is  shown  In 

Figure  3.  A  linear  variation  of  In  (<x/Kio)  indicates  a  constant  growth  rate 

for  exponential  variations  In  Kj.  The  data  are  fit  by  a  linear  least  squares 

technique  to  determine  the  growth  rate.  Three  cylinders  were  tested.  Height 
(2c),  diameter  (2a),  aspect  ratio  (c/a),  and  offset  (h)  of  the  geometric 
center  of  the  cylinder  from  the  glmbal  pivots  are  listed  below: 
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Cylinder 

Type 

Height 

(cm) 

Diameter 

(cm)  c/a 

Offset 

(cm) 

1/1 

13.246 

12,718  1.042 

-1.530 

3/2 

18.854 

12.690  1.486 

+5.887 

3/1 

25.768 

8.242  3.126 

-0.  959 

The  cylinders  are  referenced  by  their  approximate  aspect  ratio  which 
corresponds  to  the  cylinder  type  as  listed  above.  Silicon  oils  were  used  as 
test  fluids.  The  kinematic  properties  of  these  oils  are  given  below. 

Viscosity  Ratio 

( v/  Hvater^ 

Kinematic  Viscosity  (cs) 

Density  (qm/cc) 

100 

102.4 

0.972 

500 

528,6 

0.972 

1,000 

1,004.5 

0.9/4 

10,000 

9,468.5 

0.953 

60,000 

58.955 

0.960 

100,000 

101,772 

0.963 

Watar  has  a  kinematic  viscosity  of  approximately  1  cs  =  0.01  cm2/s  under 
standard  conditions.  Most  of  the  data  were  taken  at  a  spin  rate  of  50  Hz; 
however,  some  data  trials  were  made  at  33.3  Hz,  41,6  Hz,  or  58.3  Hz.  For  a 
particular  cylinder,  viscosity,  and  spin,  the  stability  of  the  gyroscope  was 
measured  over  a  wide  range  of  coning  frequencies.  For  the  higher  viscosity 
oils,  a  rapid  survey  at  three  coning  frequencies  was  used,  while  for  the  lower 
viscosity  oils,  a  dense  survey  of  coning  frequencies  was  made  to  produce  a 
detailed  response  curve.  In  summary,  the  parameters  that  were  varied  during 
the  experiments  were;  cylinder  aspect  ratio  (c/a),  kinematic  viscosity  (v), 

spin  rate  ( b) ,  and  coning  frequency  (t).  Control  of  these  parameters  produced 
the  following  ranges:  5  <  Re  <  12,000  and  0.04  <  t  <  0.1. 


III.  LIQUID-FILLED  GYROSCOPE  EXPERIMENTS 

Two  objectives  were  outlined  for  the  gyroscope  tests:  (1)  Determine  the 
yaw  moments  produced  for  low  Reynolds  numbers  and  (2)  Determine  the  validity 
of  linear  models  that  utilize  boundary- layer  corrections  (References  4  and  5) 
as  a  function  of  Reynolds  number.  The  liquid  moment  will  have  a  complicated 
behavior,  but  the  primary  independent  variables  are  Re,  c/a,  r,  Kj,  and  e. 

Previously  during  experimental  investigations  of  high  Re  rotating  flows, 
nonlinear  or  nonstationary  behavior  has  been  observed  as  a  function  of 
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Kr  1,12  In  the  present  experiments  for  Re  <  12,000,  such  nonlinear  or  non¬ 
stationary  behavior  was  not  observed  for  <  5  deg;  i.e.,  for  exponential 
behavior,  constant  growth  rates  were  observed.  The  basic  measurement  of  a  yaw 
growth  rate  will  normally  be  converted  Into  CLSM.  The  use  of  CLSM  will  also 

allow  for  the  comparison  of  the  gyroscope  data  with  existing  flight  data. 

^LSM  *s  usually  presented  as  a  function  of  the  coning  frequency  (x)  or  the 
product  Re  *  t,  which  is  a  Reynolds  number  whose  characteristic  frequency  is 
based  upon  the  fast  yaw  frequency  ($})  rather  than  the  spin  rate  ($). 

An  estimation  of  experimental  errors  for  the  gyroscope  data  was  not 
specifically  performed.  Rather,  multiple  runs  for  the  same  set  of  test  condi¬ 
tions  were  recorded  and  analyzed.  All  of  these  runs  are  provided,  thus  Indi¬ 
cating  the  repeatability  of  the  entire  measurement  and  data  reduction  proce¬ 
dure. 


The  wide  range  In  Reynolds  number  Is  primarily  achieved  by  the  choice  of 
the  viscosity  of  the  liquid.  For  Re  >  5,000,  the  aspect  ratios  of  the  1/1 
and  3/1  cylinders  were  selected  such  that  t  was  sufficiently  close  to  an 
eigenfrequency  to  produce  unstable  yawing  motion.  (The  aspect  ratio  of  the 
3/2  cylinder  was  specifically  chosen  not  to  be  close  to  an  eigenfrequency  and, 
therefore,  would  have  a  stable  yawing  motion  for  large  Re.)  Comparisons 
between  experimental  data  and  models  that  use  boundary- layer  viscous-type 
corrections  should  be  good  for  "large"  Reynolds  numbers.  Figure  4a  provides 
such  a  comparison  for  the  1/1  cylinder  when  Re  *  12,400.  Within  Reference  4, 
growth  rate  predictions  were  made  for  Incremental  changes  In  c/a.  A  slightly 
better  correlation  with  the  experimental  data  was  obtained  when  the  computa¬ 
tional  aspect  ratio  was  Increased  by  0.5%.  Figure  4a  does  not  Include  pre¬ 
dictions  of  growth  rate  for  slightly  modified  c/a  values.  Figures  4b  and  4c 
show  comparisons  of  theory  and  experiment  for  c/a  =  1/1  for  Re  =  2,400  and  Re 
=  1,260.  It  is  not  clear  whether  the  differences  between  the  observed  and 
predicted  growth  rates  are  due  to  small  measurement  errors  in  c/a  or  Increased 
viscous  forces.  However,  the  agreement  Is  still  sufficiently  good  for  practi¬ 
cal  applications.  Figures  5a  and  5b  present  comparisons  for  the  3/1  cylinder 
when  Re  =  5,210  and  Re  =  1,010.  The  agreement  between  experiment  and  theory 
is  not  good.  Within  Reference  4,  comparisons  of  theory  and  experiment  were 
substantially  better  for  Re  =  9,000  and  Re  =  520,000.  All  of  the  data  In 
Figures  4-5  have  a  similar  character:  nonmonotonic  dependence  of  growth  rate 
with  coning  frequency.  For  Re  <  1,000,  the  nature  of  the  response  curve  (for 
the  same  range  In  coning  frequency)  changes  dramatically.  Experimental  data 
for  these  cases  will  be  presented  In  the  form  of  a  liquid  side  moment  coeffi¬ 
cient  (CLJiM)  rather  than  growth  rate  to  allow  a  comparison  with  flight  data. 


11,  w.E.  Scott  and  W.P.  D'Amico,  "Amplitude-Dependent  Behavior  of  a  Liquid- 
Filled  Gyroscope, "  Journal  of  Fluid  Mechanics,  Vol.  60,  Part  4,  1973, 
pp.  751-758. 

12.  Richard  D.  Whiting ,  "An  Experimental  Study  of  Forced  Asymmetric  Ost illa¬ 
tions  in  c  Rotating  Liquid-Filled  Cylinder, "  BRL  Technical  Report  Go. 
ARBRL-TR-02376,  October  1981  (AD  A107948). 
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A.  Data  for  a  Viscosity  Ratio  of  10,000 

Figures  6a-c  give  versus  x  for  a  viscosity  ratio  of  10,000.  It  is 

seen  that  C^M  is  nearly  linear  with  x  for  all  cases.  Small  changes  in  Re 

were  obtained  by  different  spin  rates.  Changes  in  C|_SM  that  were  produced  by 

variations  in  spin  can  be  correlated  by  using  Re  •  x  rather  than  simply  t,  as 
seen  in  Figures  6d-f.  However,  CL$M  data  for  different  aspect  ratios  (but 

identical  spin  rates)  are  not  well  correlated  by  Re  •  x. 

B.  Data  for  a  Viscosity  Ratio  of  60,000 

Figures  7a-c  show  C|_SM  versus  x  for  a  viscosity  ratio  of  60,000.  As 

before,  it  is  observed  that  C^  is  linear  with  x  and  that  small  variations  in 

spin  can  be  correlated  by  using  Re  •  x  (Figures  7d-f).  However,  in  some  cases 
for  lower  spin  and/or  coning  frequencies,  the  gyroscope  was  stable.  Figure  7g 
also  shows  that  data  from  different  aspect  ratios  are  not  well  correlated  by 
Re  •  x. 

C.  Data  for  a  Viscosity  Ratio  of  100,000 

Figure  8  gives  data  for  a  viscosity  ratio  of  100,000.  Very  little  data 
were  taken  since  the  gyroscope  was  often  stable.  A  linear  variation  of  C. 
with  x  is  observed. 

0.  Correlation  of  Data  for  Viscosity  Ratios  of  10,000,  60,000  and  100,000 

A  large  volume  of  data  now  exists  for  low  Re  where  CLSM  varies  linearly 

with  x.  All  of  the  data  are  shown  in  Figure  9  with  Re  •  x  as  the  correlation 
parameter.  From  previous  plots,  it  was  expected  that  a  good  correlation  would 
not  occur.  However,  the  expanded  scale  artificially  enhanced  the  correlation, 
except  new  the  data  are  clustered  based  upon  the  viscosity  ratio  of  the  test 
liquid;  l.e.,  all  the  data  for  Re  •  x  >  4  are  for  a  viscosity  ratio  of  10,000. 
The  CLSM  data  were  plotted  against  combinations  of  (Re,  x,  c/a).  A  corre¬ 
lation  for  C^sm  3  0.0125  ♦  0.0756  (Re)( x2)(c/a)_l  was  selected  by  a  linear 

least  squares  fit  and  is  presented  in  Figure  10.  Several  deficiencies  are  ob¬ 
vious.  First,  much  of  the  data  for  very  low  Re  are  collapsed  near  the  ori¬ 
gin.  However,  the  larger  (and  potentially  destabilizing)  values  of  CL$M  are 
well  correlated. 

E.  Comparison  With  Flight  Data 

Flight  data  from  projectiles  are  available  within  References  8  and  9. 
Unfortunately,  flight  data  and  gyroscope  data  are  not  at  the  same  aspect  ratio 
and  Reynolds  number,  although  new  tests  are  planned  to  match  flight  and  gyro¬ 
scope  conditions.  Yawsonde  data  were  processed  for  C^SM  in  a  fashion  similar 

to  that  of  the  gyroscope  data.  The  amplitude  of  the  yaw  was  greater  than  that 
of  the  gyroscope.  The  effects  of  aerodynamic  damping  were  neglected,  but  the 
flight  Mach  numbers  were  transonic  and  under  such  conditions  the  aerodynamic 
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damping  of  the  projectile  is  small.  Figure  11  shows  a  comparison  between 
flight  data  (Reference  9,  Round  El-9542)  and  gyroscope  data  for  comparable 
Reynolds  numbers.  Figure  12  shows  a  second  comparison  between  flight  data 
(Reference  8,  Round  El-9396)  and  gyroscope  data.  Due  to  the  differences  in 
both  Re  and  c/a,  comparisons  may  not  be  valid.  However,  the  flight  and 
gyroscope  determined  CLSM  data  are  similar  in  both  cases.  The  data  in  Figure 

11  show  a  consistent  trend  in  t,  but  this  did  not  occur  in  Figure  12.  A  data 
comparison  for  a  wide  range  of  t  would  be  difficult,  since  projectiles  nor¬ 
mally  have  t  ~  0.08.  The  flight  data  can  also  be  compared  with  the  correla¬ 
tion  suggested  in  Figure  10.  For  Re  =  13.1  and  =  0.0153,  (Re)(T)(c/a)_1 

=  0.024,  while  for  Re  =  336  and  =  0.020,  (Re)(t)(c/a)-1  =  0.76.  The  Re  = 

13.1  datum  falls  close  to  the  suggested  correlation,  but  the  Re  =  336  datum 
does  not  agree  with  the  correlation.  The  poor  correlation  of  the  Re  =  336 
datum  could  be  attributed  to  the  differences  seen  in  Figure  12  for  the  simple 
comparison  of  versus  t. 


IV.  DISCUSSION 

Yaw  moments  from  the  gyroscope  can  be  readily  displayed  as^simple  func¬ 
tions  of  the  Reynolds  number.  For  data  sets  where  c/a,  v,  and  are  con¬ 
stants,  C[_SM  can  be  averaged  over  t  or 

CL$M  3  f  tc/a»  Re)  * 
ave 

Figures  13a-c  show  C.  <vM  versus  Reynolds  number.  The  scales  of  these  fig- 

Lonave 

ures  are  Identical  to  emphasize  various  aspects  of  the  data.  In  Figure  13a, 

Ci  cu  dramatically  Increases  as  the  Reynolds  number  increases.  C,  cM  does 
Lonave  ave 

not  decrease  for  extremely  high  Re  values  (as  did  the  despin  moment  data  in 
Figure  1).  The  resonance  between  the  coning  frequency  and  the  liquid  eigen- 
frequency  for  c/a  =  1/1  only  becomes  stronger  as  Re  increases.  A  similar 
situation  exists  for  c/a  =  3/1,  However,  the  effects  are  much  less  obvious, 
as  seen  in  Figure  13b.  The  resonance  with  the  liquid  eigenfrequency  is  much 
weaker  for  c/a  »  3/1  than  for  c/a  =  1/1.  The  3/2  cylinder  was  specifically 
chosen  so  as  not  to  produce  a  resonance  within  the  range  of  experimentally 
produced  coning  frequencies.  The  yaw  behavior  for  c/a  =  3/2  Indicates  a 
maximum  C.  cM  for  Re  <  100,  with  smaller  yaw  moments  for  either  larger  or 

i-^ave 

smaller  Reynolds  numbers  (as  in  Figure  13c).  This  type  of  behavior  is  similar 
in  character  to  the  despin  moment  data, 

A  typical  response  curve  for  large  Reynolds  numbers  is  shown  in  Figure 
4a.  Response  curves  for  lower  Reynolds  numbers  show  substantially  reduced 
growth  rates  (Figures  4b  and  4c).  The  predicted  growth  rates  increase  linear¬ 
ly  with  t  until  the  maximum  growth  rate  is  achieved.  Also  note  from  Figures 
14a-c,  that  as  the  Reynolds  number  was  reduced,  the  coning  frequency  for 
maximum  growth  rate  increased.  For  Re  <  1,000,  the  maximum  growth  rate  may 
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shift  to  t  >  0.1.  Thus,  growth  rates  would  increase  linearly  with  coning 
frequency  for  t  <  0.1.  This  is  exactly  what  is  observed.  A  continued  reduc¬ 
tion  in  Re  will  eventually  give  extremely  large  viscous  damping  and  unstable 
yaw  behavior  will  not  be  possible. 

The  previous  rationale  is  not  applicable  to  the  3/2  cylinder.  Murphy4 
has  shown  that  the  average  value  of  the  liquid  moment  is  increased  by  a  reduc¬ 
tion  in  Re.  Hence,  for  a  nonresonant  cylinder,  such  as  c/a  *  3/2,  large  Re 
behavior  is  stable  while  low  Re  behavior  could  be  unstable.  When  Re  approaches 
unity,  however,  the  liquid  will  act  as  a  rigid  body  and  no  destabilizing 
effects  will  be  observed.  This  mechanism  is  in  qualitative  agreement  with  the 
data  for  the  3/2  cylinder  as  shown  in  Figure  13c. 


V.  SUMMARY 

Comparisons  between  gyroscope  data  and  theory  were  consistent  for  c/a  = 
1/1  for  Reynolds  numbers  as  small  as  1,010,  but  comparisons  were  Inconsistent 
for  c/a  »  3/1  when  Re  *  5,210.  This  Indicates  that  theoretical  models  with 

-1/2 

viscous  corrections  of  order  Re  m^y  have  a  Reynolds  number  and  aspect 
ratio  limit  when  the  Reynolds  number  approaches  1,000.  Yaw  moments  for  highly 
viscous  flows  are  probably  not  generated  from  a  new  mechanism,  but  rather  are 
the  remnants  of  high  Reynolds  number  resonances.  However,  the  liquid-induced 
yaw  moments  are  sufficiently  large  for  10  <  Re  <  100  to  destabilize  the  yawing 
motion  of  a  spinning  projectile  or  a  gyroscope. 


typical  Response  Function  for  a  High 
teynolds  Kurrfcer,  Liquid-Filled  System 
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Comparison  of  Theory  by  Murphy  (Ref.  4)  and  Gyroscope  Eafca:  Fill  Ratio 
100Ss  c/a  «  3.126,  Re  =  1,010,  I  =  10.76x10s  gro^cm2,  n  «  ls  j  =  1 


Figure  6 d.  liquid  Side  Moment  Coefficient  (CLSM)  Versus  Re  •  t  for  c/a  of  1/1 
Kith  a  Viscosity  Ratio  of  10,000 
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Reynolds  Number  •  Con i ng  Frequency 

Figure  6e.  Liquid  Side  Moment  Coefficient  (C^gjjj)  Versus  Re  »  t  for  c/a  of 
3/2  Mith  a  Viscosity  Ratio  of  10,000 
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Figure  7a.  Liquid  Side  Moment  Coefficient  (C^)  Versus  Coning  Frequency  {%) 
for  c/a  of  1/1  and  a  Viscosity  Ratio  of  60,000 
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030r  Symbol  c/a  Re  Spin  (Hz) 
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Reynolds  Number  •  Coning  Frequency 

Figure  7d.  Liquid  Side  Moment  Coefficient  (C|_SM>  Versus  Coning  Frequency  (t) 
for  c/a  of  1/1  and  a  Viscosity  Ratio  of  60,000 
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Reynolds  Number  •  Coning  Frequency 

Figure  7e.  Liquid  Side  Moment  Coefficient  (CLSH)  Versus  Re  •  x  for  c/a  of  3/2 
and  a  Viscosity  Ratio  of  60,000 
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Figure  7f.  Liquid  Side  Moment  Coefficient  (CLSM)  Versus  Re  * 
for  c/a  of  3/1  and  a  Viscosity  Ratio  of  60,000 
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Figure  7g.  Liquid  Side  Moment  Coefficient  (CLSM)  Versus  Re  •  x 

for  c/a‘s  of  1/1,  3/2,  and  3/1  With  a  Viscosity  Ratio  of  60,000 
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Correlation  of  all  Data  with  Viscosity  Ratios  of  100,000,  60,000,  and 

10 ,000  for  C.  Versus  Ee  *  t 
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LIST  OF  SYMBOLS 


radius  of  a  right-circular  cylindrical  cavity  containing  liquid 

liquid  in-plane  moment  coefficient  for  one-mode  coning  or  spiral 
motion;  the  imaginary  part  of  Cj_M 

fast  (j=l)  and  slow  (j=2)  mode  liquid  in-plane  moment  coeffi¬ 
cients;  the  imaginary  part  of  Cyj. 


fast  (j=l)  and  slow  (j=2)  mode  liquid  moment  coefficients 


mLa2<}>2T1K12CLRM 


liquid  side  moment  coefficient  for  one-mode  coning  or  spiral 
motion;  the  real  part  of 

fast  (j=l)  or  slow  (j=2)  mode  liquid  side  moment  coefficients; 
the  real  part  of  C^. 


I  Magnus  moment  1 
(1/2)  PSt2  V  i  U| 

[sum  of  the  damping  moments 
(1/2)  pSt 2  V  |cross  spin) 


I  static  moment! 
(1/2)  pSt2  V  \l\ 


(normal  force) 
(1/2)  PSV2  l 


one-half  the  length  of  the  cylindrical  cavity  containing  liquid 


(pS&/2m)  (V/A) 


CN  -  kw"2  f CM  +  C, 


(vCmJ 


a  y  T  M. 
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distance  from  the  gyroscope’s  pivot  location  to  the  center  of 
the  cylindrical  cavity 

axial  and  transverse  moments  of  inertia  of  the  gyroscope  or 
projectile 

magnitude  of  the  j-th  yaw  arm  (j=l,  2) 
initial  value  of  Kj 

(Ix/m£2) 1/2,  the  projectile's  axial  radius  of  gyration 

(Iy/mi2)1/2,  the  projectile's  transverse  radius  of  gyration 

reference  length 
(pS£3/2i  )  (v/£)2  CM 

components  of  the  aerodynamic  moment 
projectile  mass 

liquid  mass  in  a  fully  filled  cylindrical  cavity 

a2*$/\>»  Reynolds  number 
reference  area 

o24>2/ 4M,  the  gyroscope  stability  factor 
(PS£/2m)  (V/£) 

time 

magnitude  of  the  projectile's  velocity  vector 
nondlmensionallzed  growth  rate  for  single  mode  motion,  j-1 


\ *  k*'2  \ 
m  ^  , 
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ej  nondlntensionallzed  growth  rate  of  the  j-th  yaw  mode  (j*l,  2) 

v  kinematic  viscosity 

•N* 

5  complex  yaw  in  the  nonrolling  coordinate  system 

p  air  density 

0  Ix/Iy 

T  tj  for  one-mode  yawing  motion 

Tj  the  nondimensionalized  frequency  of  the  j-th  yaw  mode 

(j«l»  2) 

♦  *  t 

<t>j  (j>j0  +  <}i  t  (j*l,  2) 

tjo  initial  orientation  angle  of  the  j-th  yaw  ara  (j=l,  2) 

t 

$  spin  rate  (assumed  to  be  positive) 


Superscripts 

(♦)  vector  quantity 

(*)  time  derivative 
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Mr.  McKivrigan 
Nr.  F.  Dagostln 
Mr.  C.  Hughes 

Director,  USACSL 
Bldg.  E3516,  EA 
ATTN:  DRDAR-CLB-PA 

Mr.  M.  Hiller 
DRDAR-CLB-PA 
DRDAS-CLN 
DRDAR-CLJ-L 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  beiow;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1.  BRL  Report  Number _ 


2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  <r  procedure,  management  procedure,  source  of 
ideas,  etc.) 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  oe  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 


Name: 

Telephone  Number: 
Organization  Address: 


